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Of the two main objectives we follow in this paper, the first one consists in the studying operators of the 
form (3, - iAr)“, a = 1/2, -1/2, -1,..., where Ap is the Laplace-Beltrami operator. These operators arise in 
the context of nonreflecting boundary conditions in the pseudo-differential approach for the general Schrodinger 
equation. The definition of such operators is discussed in various settings and a formulation in terms of fractional 
operators is provided. The second objective consists in deriving comer conditions for rectangular domains which 
are fundamentally inadmissible in the pseudo-differential approach. Stability and uniqueness of the solution is 
investigated for each of these novel boundary conditions. 


I. INTRODUCTION 

In this article we consider the problem of construction of 
nonreflecting boundary condition for the general two dimen¬ 
sional Schrodinger equation. In particular, we consider the 
following initial-value problem (IVP): 

idfU + Au + (p(x, t, \uf')u — 0, (x, f) e X R+, 
m(x, 0) = Mo(x), X 6 R^. 

The initial data is assumed to be supported within the com¬ 
putational domain, Q,-, i.e., supp^Mo(^ c Q,-. This problem 
has been treated by several authors lUl- U^ . Exact formula¬ 
tions of the transparent boundary condition (TBC) for the free 
Schrodinger equation on convex domains, Q,, with smooth 
boundary was provided by Schadle [@1 in terms of a single and 
a double layer potential. The special case of a circular domain 
was treated by Han and Huang Js]]. Earlier attempts to derive 
an exact TBC for a rectangular domain by Menza proved to 
be problematic on account of the presence of corners cm. 
This problem is resolved in a recent work by Eeshchenko and 
Popov il. On account of the lack of integrability, these tech¬ 
niques cannot be applied to the general Schrodinger equation 
and one has to turn to approximate methods (refer to the re¬ 
view article Cl for a comprehensive literature survey). 

Eor the approximate methods, we restrict ourselves to the 
pseudo-differential approach (in particular, the gauge trans¬ 
formation strategy) for constructing approximate nonreflect¬ 
ing boundary conditions referred to as the absorbing bound¬ 
ary conditions or artificial boundary conditions for various 
types of computational domains. Our goal is to understand 
operators of the form idt - iA^Y, a - 1/2--1/2,-1,..., 
which a|meared in the works of Menza Cl H and Antoine 
et al. Sd- Several aspects of such artificial boundary con¬ 
ditions (ABCs) which comprises these operators are not quite 
well understood; we discuss these issues which motivate the 
present work in the subsequent paragraphs. 

Contrary to the existing belief that the aforementioned op¬ 
erators can only be implemented via a Pade approximation 
of a monomial of fractional degree a, i.e., z“, it is shown in 
various settings that this operator can be expressed in terms 
of fractional operators. Eor arbitrary functions f{x, t), {x, t) 6 
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RxR+, the operation ( 5 ,-/Aj:)“/(a, f). a — l/2,-l/2,-l,..., 

requires the knowledge of the function over its entire support. 
If supp^/(x, t) c Ex, Vf > 0 or f{x, t) is periodic with respect 
to X and it can be uniquely defined by its values at .r 6 r;^ for 
allf > 0, then it is shown that a formulation particularly con- 
venient for expressing the TBCs/ABCs for the IVP in ([TJ can 
be developed^ Let us remark that the numerical implemen¬ 
tation of such operators is not being considered in this paper; 
however, it can be easily seen that the formulation developed 
in this paper makes these operators amenable to convolution 
quadrature. It might be expected that the new scheme affords 
improvement in accuracy over the existing Pade approxima¬ 
tion based method0 reported in ina. 

Eurther, it is well known that the pseudo-differential ap¬ 
proach cannot be applied to computational domains with cor¬ 
ners. This precludes the rectangular domain which happens 
to be a very convenient choice of the computational domain. 
The ABCs involving the operators of the form (5, - /Ap)" can¬ 
not be adapted to the rectangular domain in a straightforward 
manner; however, the ABCs obtained as local approximations 
(with respect to x) or equivalently the high-frequency approx¬ 
imations of this operator admit of the possibility of construct¬ 
ing corner conditions. We demonstrate this possibility for 
the free as well as the general Schrodinger equation given by 
(fTb. Our approach is closely related to the ideas presented 

inlUlIl. 

Another program that we have followed in this paper is of 
obtaining the well-known energy estimate 

I|m(x, 011^(0;) < l|Mo(x)||L2(n,), (2) 

for the IVP in ([T]| with TBCs/ABCs as boundary conditions 
involving operators of the form (3, - iA^Y or high-frequency 
approximations of it0. Under the assumption that the solution 
exists, this inequality guarantees the stability as well as the 


* The TBCs derived by Feshchenko and Popov (3l happen to be a special 
case where this operation can be carried without having to consider the 
function f{x, t) over its entire support. The function f{x, t) here refers to 
the restriction of the field w(x, t), x € to the segments of the boundaiy 
of the rectangular domain. 

^ A distinction must be made between the Pade approximation based meth¬ 
ods which are applied directly to the operator {dt - those applied 

to the fractional operator formulation of this operator discussed in this pa¬ 
per. 

^ The corner conditions are also incorporated in this equivalent formulation. 
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uniqueness of the solution of the equivalent initial boundary- 
value problem (IBVP). In certain cases, this result can be ob¬ 
tained by resorting to a general construct such that for any 
pseudo-differential operator, P, and a function u 6 C^(f2) we 
have 

2 ‘^{u\Pu) - {u\Pu) + {u\P^u), ( 3 ) 

where is adjoint of the operator P, stands for real part 
and {u\v) - u*vdQ. Given the symbol crp of P, the symbol 
of the adjoint, cr^t, can be computed using the following gen¬ 
eral formula: Assuming y e R" and ^ the covariable of y, then 

crMy, 0 ~ X 

aGN" 

where cr^ stands for complex conjugate of crp (complex con¬ 
jugate of z 6 C is also denoted by z). A more detailed discus¬ 
sion of this approach is provided inlAl 

The discussion of the primary results in this paper is 
broadly divided into two sections: Section |II] deals with the 
free Schrodinger equation while Section |III] deals with the 
general Schrodinger equation. For each of these problems, 
we consider two types of domains, namely, the rectangular do¬ 
main (or, infinite strip with periodic boundary condition along 
the unbounded direction) and convex domains with smooth 
boundary. The basic definition of the operator (dt-iAp)", a - 
1 / 2 , - 1 / 2 , - 1 ,..., is discussed in various settings. Two fam¬ 
ilies of ABCs are considered in this paper: first one obtained 
via the standard pseudo-differential approach and the second 
one obtained as the high-frequency approximation of the for¬ 
mer. The derivation of corner conditions and the study of sta¬ 
bility and uniqueness of the solution are carried out separately 
for each of these problems in the subsections. 

II. FREE SCHRODINGER EQUATION 

Let us start our discussion with the linear case of the equa¬ 
tion ([T]i with null-potential, i.e., 

id,u + Au - 0, (x, f) 6 X R+, 

( 5 ) 

m(x, 0) = Mo(x), X 6 R^. 

The problem of constructing transparent boundary condition 
can be exactly treated for the case of compactly supported ini¬ 
tial data for a computational domain Q, with a general smooth 
boundary T i]. The basic approach involves solving the ini¬ 
tial boundary value problem on the exterior domain. The ex¬ 
terior domain is defined by Qf = R^ \ Q,. It is known that a 
boundedness condition at infinity does not ensure unique so¬ 
lution of the above IVP; one additionally needs to impose a 
Sommerfeld-like radiation condition at infinity in order to ex¬ 
clude all the incoming waves from infinity. This condition 
reads as 

■ -^1 + e^‘*dy^uj- 0. (6) 


Using the decomposition of m(x, f) 6 L^(R^) = L^(G,)©L^(f 2 f) 
and introducing the fields v(x, t) and w(x, t) we have 

idfV -H Av = 0 , (x, t) e £ 2 / x R+, 
v(x, 0 ) = Mo(x), X 6 G;; 

id/W + Aw = 0 , (x, f) 6 Ge X R+, 
u(x, 0 ) = 0 , X e Q(., 

limixHoo VW^Aw ■ || -I- e’'?5y^wj= 0; 

v(x, Olr = vv(x, f)lr, d„v(x, f)lr = 5„w(x, f)lr- 

Construction of the nonreflecting boundary conditions in¬ 
volves solving the exterior problem exactly. In the follow¬ 
ing sections we first consider a infinite strip as computational 
domain (or periodic boundary condition in the direction it ex¬ 
tends to infinity) then extend the results to a rectangular do¬ 
main. 
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FIG. 1. The computational domain in the form of an infinite strip 
is depicted here. This domain can be replaced by a rectangular do¬ 
main if we assume periodic boundary condition along the unbounded 
direction. 


A. Infinite strip 

In this section we restrict ourself to the case of an infi¬ 
nite strip with the boundary parallel to one of the coordinate 
axes, or, rectangular domain with a periodic boundary con¬ 
dition along one of the coordinate axes. Let the coordinate 
axes be labelled as Ci and 62. For the infinite strip, say with 
boundary parallel to the axis 62, we assume that supp^.^ uo(x) is 
bounded. The derivation of the TBCs for the infinite straight 
boundary is particularly simple and can be obtained using 
Laplace transform in time and a Fourier transform in space. 
Let us consider the IVP in Q where the domain is defined 
by the infinite strip between xi - x; and xi = x^. For the 
TBC on the right boundary, one must consider the IBVP on 
the right exterior domain Q; = [x^, cm) x R. Let the covari¬ 
able of (xi, X2) be denoted by (/"i, (2), respectively. We intro¬ 
duce the notation ^xif{x\,X2,t) - ■^xi{f}{^\,X2,t) for one 


lim 

IxHoo 
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dimensional Fourier transform with respect to xi (similarly 

^ 2 j 0 for Fourier transform with 
respect to xi). For denoting the Laplace transform of a func¬ 
tion we use = ^,[f]{s) = F{s). 

Let us denote the Fourier transform with respect X 2 of 
w(xi, X 2 , t) by vv(xi, ^ 2 , 0 and Laplace transform with respect 
to t of vv(xi, X 2 , f) by ^(-* 1 : ^ 2 , ■s), i-e-, 


where represents the convolution operation, we have 




id,w(x\_,i^ 2 , t) + {dl w){x]_,^ 2 , t) 


( 10 ) 


W(X1,^2, i) = ^x2WiXx,X2,t), 

( 8 ) 

W{X\,^2,S) - ^tW(x\,^2,t)- 

It is also easy to verify 

For the case of compactly supported initial data, we have 


+ a^)W{x\,^ 2 ,s) = 0, XI 6 [Xr, oo), (9) 

_ dte~'^^*d^^^^e'^^'w{x\,^2,t) - e~‘^^'d^^^^e'^^'dtw(xi,^2,t)- (H) 

where a - Jis - such that 3 (a) > 0. The solution can be 
worked out as follows: observing 

transform, we obtain the 


dx,W(x\_,(2,s) = /q'VT(xi,^ 2, s), 

jSf-'[5„lT(xi,^2, s)] = * J^'‘[a^lL(xi,^2, «)], 


Now, taking the inverse Fourier 
Dirichlet to Neumann map as ^ 


dx,w{s.,t) = 


„ot/4 nt n 

V)r Jo Jr 


idjwixi , xL r) -H d\ w(xi, xL t) 


^(X2 - xl, f - t) 


VT 


dx'^dT 


= -(5, - iSb- 


g-/;r/4 nt n 

Vtt Jo Jr 


&(x2 - xLt - t) ^ 
w(xi, x^, t)- -=^ - dx2dT, 

xjt - T 


( 12 ) 


where the convolution kernel is given by 

p—injA ^ 

G{x 2 , f) = , with #(^ 2 ,0 = (13) 

V47rf 

This map can be expressed compactly if we employ the nota¬ 
tion 


Next, we would like to obtain a more local approximation 
of this boundary condition valid for small times. To this end, 
let us consider 





(17) 


r r f(x\,x 2 ,T) 

yTT Jo Jr 


^(X2 - x'2,t- r) 

xjt - T 


c/x^t/r, 
(14) 


where / is assumed to be an element of C“(R x R+) such that 
it is of the Schwartz class with respect to X 2 . On account of 
the singular nature of the symbol, it is not a pseudo-differential 
operator. Flowever, away from the points satisfying = 0, 
it can be microlocally regarded as a pseudo-differential oper¬ 
ator with a symbol (i^ + with the branch cut defined 

by -7T < arg{i^ + ^ where {(,2,0 the covariables 

of (x 2 , t). Further, if we take the operator (J, - to be 

defined by 


Setting ( - St, we have 

— r 

TTsJi Ja¬ 


in Vf 
1 




--d( 


(hym 


r j-(i _ hi] 

Ja-icc V^v t 


In s/i Ja-, 

1 r"'“/ 1 

rV^ «Jfl—I'oo \ V? 


dt 


2n s/i 
yK/4 


4- 


e^d( 


(18) 


+ 


2(i^)3/2 8(if)5/2 

36^/4 ,,2 . 

— r/ t2 " 


... y^d( 


mi)' 


sni) 


{d, - idif^f = {d, - idim, - idiy/^f], (15) 


then 

5x,w(x,f)+e""'^'‘(5,-i5j^)'''V(x,f) = 0, (x,f)eO;xR+. (16) 


Taking the inverse Fourier transform, we obtain the asymp¬ 
totic expansion as 
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For the periodic case, we may take the computational do- It is also easy to verify the following relations: 
main to be Q,- = [x/, x^] x [ 0 , Itt] so that 


w(xi,X 2 ,f) = ^ W;„(xi,f)e™''J 

(20) 



meZ 




W(xi,X2,i) = ^W;„(xi,i)e™"^ 

(21) 



meZ 


^ e-‘'^'‘d;^^^f'^"'d,Wm(xi,t), 

(23) 


with Wm{x\,s) - ^[vvm(xi,f)]- For simplicity we demon¬ 
strate the procedure for the m-th Fourier component. The 
complete result is obtained by superposing all the compo¬ 
nents. Putting am = Vii - the following results can be 
obtained in the same manner as before to obtain 


-inrt 


id,Wm(xut) - irrwmixiJ) 


^ ie-''”"‘dy^e'"'"‘wm(xut). 


(24) 




id,Wm(xi,t) - rrrwmix\,t) 


. ( 22 ) 


Using these relations, one can derive an equation analogous to 

(Ha 


^/;r/4 nt 

5j,w(x, f)= —— I I /( 9 tW(Xi,X 2 ,t) -H ( 9 y w(xi,X 2 ,t) 
Jo Jo L 2 


g{x2 - xi, f - t) 




t/xjc/r 


- -{d, - idi)- 


g-OT/4 

Vtt Jo Jo 


Q{x2 — x^, f — t) ^ 

w(xi, Xj, t)- —= - dx2dT, 

yjt — T 


(25) 


with the convolution kernel given by 

g{x2,t)^Yj^’"""\ (26) 

meZ 

which is defined only in a distributional sense. For periodic 
functions the operator defined in (fl4l l becomes 


f f f(Xl,x'2,T) 

yn Jo Jo 


g(x2 - xl, f - t) 


VT 


t/X2<iT, 

(27) 


where the kernel is defined by (l26T l. We now turn our attention 
to obtaining a form of the operator (5, - which can be 

numerically implemented. Using the relation (l24l i. we may 
write 






1/2 


^w„(xi,f)e" 




meZ 


(28) 




Introducing the auxiliary function ip(xi,X 2 , t, t') such that 


^(xi,X2,f,f') = J][e-™'''-''V^(xi,f')J”’"^ 

meZ 

J^,w(x, t) = -^2, t, f')|r=,. 


(29) 


In order to determine all the values of the function needed 
to compute the non-local fractional derivative, the following 


relations can be used: 

[idt + d\^]ip{x\, X 2 , t, f ) = 0, xi e [xr, oo), 
ip(Xi,X2, f , f) = w(Xi, X 2 , t'). 


(30) 


For all possible values of f satisfying 0 < f < t and X 2 6 
[0,27r], one must solve the above initial value problem starting 
from the initial condition at f = f' under the periodic boundary 
condition in X 2 . This process is schematically shown in Fig. 
|2] The history of the held is needed along the vertical line up 
to the diagonal in the (f, f')-plane. 

Remark 1. Such a procedure can be used to compute the ac¬ 
tion of any operator of the form (di - i n - 1 , 2 ,..., 

on any function periodic in x e'R, say, /(x, /)■ Introducing an 
auxiliary function (p(x, t, f) such that 


^(x,f,r') = yL'™^<'-''V,„(f')J'“ , 

meZ 

(d, - idly'^^fix, t) = d^f''^ip{x, t, t')\f=,. 


(31) 


The associated IVP is given by 

[idt + i9^]i,£’(x, f, f) — 0, 
ip{x, f, f) - f(x, f). 


(32) 


1. Stability and uniqueness 

An equivalent formulation of the IVP Q on the computa¬ 
tional domain Q, = [x;, x^] x [0,27r] with periodic boundary 
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FIG. 2. Schematic illustrating how the auxiliary equation OOb will 
be solved in order to provide the history of the field needed in the 
TBCs. The field is known along the diagonal which serves as the 
initial conditions to obtain the values of the field needed in the TBCs 
(arrow in the line depicts the evolution direction in time). 
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FIG. 3. Rectangular domain 


B. Rectangular domains 


condition along the axis ea is given by 


Let us consider a rectangular domain given by 


id,u + Au -0, (x, f) 6 Q; X R+, 

m(x, 0) = Mo(x) 6 supp Mo(x) C Qj, 

u(xi, X 2 + 27T, t) — u{x\,X 2 ,f), t > 0, 

dnu(x, t) + t) -0, x 6 L/ U L^. 

Assuming that the solution m(x, t) exists, we have 

I {dt\u'^)(fx - 2% f (u*iVu) ■ dg, 

Jn, JriuFr 

I|m(x, 7’)||Ln,) - l|Mo(x)||^ 2 (n,) = ^ f f {uiVu) ■ dg, 

Jo Jr,ur, 

(34) 


Q; = {(Xi,X2) : Xi, < X2 < Xi,Xi < Xi < Xr), (38) 

and denote the boundaries as ^ = {(xi,X 2 ) e dQ., : xi - xi^r), 
respectively and L^,,, = {(xi,X 2 ) e (90,- : xi = Xhi}, respec¬ 
tively. Assuming supp,jMo(x) bounded in Q,-, the TBCs for the 
infinite strip cannot be taken to be the transparent boundary 
operators at the straight edges because the corresponding op¬ 
erator requires knowledge of the entire support of the field 
along the tangential direction at the boundary. This clearly 
cannot be achieved because after a certain time the field would 
have spread outside the domain defined by any segment of the 
rectangular domain, say, T^. This issue can be resolved in the 
following way: Using the same notation as in the last section 
and observing that 


where dg — e„|(ix| for x e T. Further noting Um{x\,t) e 
//*^^([0, T]), we have 


Ir 


--2% J" dt J" {u*.^d, - idl^u^dx 2 . 




V r dt(um(Xr,t)e ‘'"'‘*](dy^u,„(Xr,t)e 


iv r , 


dm 


1/2 




Um(Xr, t)e 


{Xr,^ 


< 0 . 


Similarly, it can be shown that 

II = -2% f dtj ^d, - idlujdx2 < 0. (36) 
Therefore, we have 


|m(x, T)||i,2(n,) < ||Mo(x)||/.2(n,). 


(37) 


This result also guarantees the uniqueness of the solution. 


djc,w(xi,mo 

d,,w(xumt) = 

djt,w(xi,X2, t) - -e~'’’^‘^d]!^if{xi , X2, t, f) 

t'=t 


(39) 


where the auxiliary function is defined by ^X 2 ‘p(xi, X 2 , t, f) = 
e~'&‘^*'^w(xi, m t')- It is easy to verify that the function ip 
satisfies the following IVP 


[id, -t- d\]ip{xi,X 2 , t, f) = 0 , 
(p(Xi,X 2 , t', t') - w( 2 i:i, 21:2, t'). 


(40) 


The boundary conditions at the endpoints of the segment F^ 
are not known because the original problem is defined on the 
infinite domain 2 C 2 6 R. However it can be shown that the 
value of the auxiliary function (f{x\,X2, t, f) at f = 0 is com¬ 
pactly supported in F, for xi e [xr, 00 ). This would allow 
one to impose the transparent boundary conditions at the end¬ 
points of F^. To this end, let us consider the IVP defined 
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by 0, it can be solved using Fourier transform in (xi,X 2 ). 
Putting Mo(^) = we have 

ip(xi,X 2 ,t,t') ^ f 

(27ry Jr2 

Therefore, supp^^i,c)(;ci, X 2 ,0, f') c F^ for xi e R and 

[idr + d^Jip(xuX2,t,t') - 0. (42) 


we obtain 


{d,-idlr"'/^f(x,t)^ 



dx'dt''K(x-x',t-t')f(x',f). (47) 


From this definition it is clear that the operator cannot be 
defined on a compact domain with respect to x for arbitrary 
functions. Let us introduce ip{x, t, t') defined by 


(f{x, t, f) - ^ JJ' dfdx'e ^ ‘ ^f(x', f), (48) 


Hence, the boundary condition at the endpoints of F^ is given 
by 

dx2‘f(xi,X2,t,f) + e~‘’'^^d]^^ip{xx,X2,t,f) - 0, (43) 

xi 6 Fr, X 2 e {xt, Xb}, respectively. This requires the knowl¬ 
edge of (p(xi,X 2 ,t,t') for X 2 - OFr and 0 < t < t'. The al¬ 
gorithm can be explained by means of the Fig. |4] There two 
IVPs for the auxiliary field, (p{xi ,X 2 ,t, t') each of which evolve 
the field either above or below the diagonal in the (f, f')-plane 
starting from the values at the diagonal in their respective do¬ 
mains. These points are represented by filled circles and the 
arrows denote the direction of evolution. All diagonal points 
are evolved along t or f'-axis in order to provide the history 
needed for TBCs in the current time step (represented by hor¬ 
izontal or vertical line in (f, f')-plane). The empty circles cor¬ 
respond to the corner points. The values of the auxiliary field 
at the corners are needed for the TBCs satisfied by the aux¬ 
iliary fields; this relationship is depicted by the broken lines. 
Note that these values at the corners can be taken from the 
adjacent segment of the boundary where it is already being 
computed. 


so that 

[d, - idlT’"^^fix, t) = 5,7'”^^v?(x, f, f')lf'=f, 

idt(p{x, t, t') + d\(p{x, t, f) - 0, (49) 

ip{x, f, f) = f(x, f). 

The possibility of defining this operator on a compact domain 
with respect to x depends on the initial condition ip{x, 0, t') 
given by 

ip(x,0,t') - — JJ' dfdx'e~‘^^^~^y'^^' fix',t'). (50) 

If suppj.ipix, 0, t') is compact then the PDF in (I491 l can be 
solved on any domain Q such that supp^^ipix, 0, t') c Vf' > 
0. This allows us to impose TBCs for (gUl at the boundaries; 
however, the value of the function tp(x, t, t') at the boundaries 
for t e [0, f'] needs to be computed. 

1. High-frequency approximation 


Remark 2. From the discussion above we obtain the follow¬ 
ing useful definition of the operator (df—id^y^^, (x, f) e RxR+.' 

idr-idiy/^fix,t) 

(44) 

This definition makes it explicit that one has to consider the 
function fix, t) over its entire support with respect to x 6 R 
in order to compute the expression; however, in special cases 
this can be avoided. Similarly, a formal definition of the op¬ 
erator idt - idy)^"‘^^,ix,t) 6 R X R+, m > Q,for any arbi¬ 
trary function fix, t) with sufficient smoothness property can 
be given. This operator can be defined as 

id,-idir'"/^fix,t) 

(45) 

The integration with respect to f can be easily performed. 
Defining the convolution kernel 'Kix, t) by 


The high-frequency approximation affords the possibility 
of simplifying the TBC by making them “local” in terms of 
the spatial variable. The asymptotic expansion worked out in 
( fT9b can be carried out for each of the edges of the rectangular 
domain to obtain the following ABCs: 

dnU + e-‘''^^dy\ - d^'^u = 0, X 6 F, U Ti, 

f (51) 

d„u + e-"l‘^d\'^u - e'”l‘^-d\d^^l^u = 0, x 6 F* U F,. 

These boundary conditions become problematic at the corners 
of the rectangular domain. This aspect can be illustrated by 
the considering the weak formulation of the original IVP as 
follows: Consider a test function fix) e //^(Q,); taking the 
inner product with the equation 0 we have 

I iidfU -H 'V^u)fd^x 

Jn, 

- id, I ufd^x- I (SIu)■ (yf)d^x-¥ I f(yu)-dq. 

Jn, Jn, Jdn, 

(52) 



.,-inlA m-l-3 

r 

‘7C(x, t) = < 

2V5Fr(,«/2)f ^ 

['47 J 


[o. 



t > 0, 
t < 0, 


Let us consider the top and right boundaries. The boundary 
(46) integrals are given by 
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t' 



t' 



FIG. 4. Schematic depiction of the evolution problem for the auxiliary field ip{xi,X 2 , t, f) in the (t, f')-plane. The filled circle corresponds to 
the entire domain while the empty circle denotes the end points. The evolution over the interior of the domain is carried either above or below 
the diagonal starting from the diagonal values (direction of evolution is depicted by the arrow head). The end points are evolved on either side 
of the diagonal. Note that the vertical/horizontal lines where the arrows end corresponds to the history of the auxiliary field needed for TBCs 
in the current time-step. The TBCs for auxiliary field require the history of the auxiliary field at the endpoints, depicted by broken lines; these 
values are taken from the adjacent segment of the boundary. 


Ir + It 



- f {dx2il^){dx2d~^^^u)dx2 + ipdx^dj 

X2=Xb JVr 


l/2„ 


Xr 

Xl-Xi 


I 


^'^u)dxx 


(53) 


The following terms corresponding to the top-right corner in 
the above equation is problematic as the boundary conditions 
in the current form cannot be used to evaluate them: 

[dx2d^^'^u +^ dJ^^^(d2C2U +djcu] .(54) 
V ''Cnr, V ''r,nr, 

In order to evaluate this term, we proceed as follows. Carrying 
out the fractional integration, of the evolution equation 
©, we have 

idy\ + {dy + dy)d;^^^u = 0, (XI,X2) 6 r, n r,. (55) 

Here, the fact that the field is zero at the corner at f = 0 is 
explicitly used to arrive at the fractional derivative. Using BCs 
(fSTb and the last equation, we obtained the following corner 
condition: 

-I- d 2 C 2 U + = 0, {xuX 2 ) 6 T, n T,. (56) 

A similar procedure can be used to the construct corner con¬ 
ditions for the other corners of the rectangular domain: 

d„u\T.+d„u\r. + -e-‘”l'^dy^u^0, (xi,.r2) e r, nr,-, (57) 

' ‘ 2 

where i + j and /, j 6 {r, f, /, b]. 


2. Stability and uniqueness 


Let us write the equivalent formulation of the IVP © for 
a rectangular domain Q,- using the TBCs derived in the last 
section: 

idtu + Au-0, (x, f) e fl,- x R+, 
m(x, 0) = Mo(x) 6 L^{D.i), suppMo(x) c Q,-, 
d„u(x, f) + e~'’'^‘^{dt - idyy^^u(x, f) = 0, x e T; U T^, ^ 
d„u(x, t) + e^'^^^'idt - id\y^^u{x, t) -0, x e T;, U T,. 
Assuming that the solution m(x, f) exits for t 6 [0, T], we have 

I|m(x, 7’)||^2(n.) - l|Mo(x)||^2(n,) 




f dt f {u*iXu)-dg+ f {u*Nu)-dg . 

Jo Ur,ur, Jrjur, 

(59) 


In the above equation, the fields in the boundary integral can 
be extended to R+ without changing the value of the integral 
by setting them identically to zero outside the interval [0, T]. 
For the right boundary, we have 

Ir = -2)R y/^t - idlu'jdx2. (60) 
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Again, without changing the value of the integral with respect 
to xi the fields can be extended to whole line by setting them 
identically to zero outside r^. We these considerations, one 
can write 


lu 


- -TK f dt f Jdt - idl^u\dx 2 , 

Jr+ Jr'' '' ' 

= J' dt J'd^ 2 (^u(xr,^ 2 ,t)e~‘^ 2 ''jl^sy^u(xr,^ 2 ,t)e^‘^^^'j, 

= r d^ f dm^ 


u{Xr,^ 2 ,t)e 


< 0 . 


(Xr,^2,0 

( 61 ) 


Similarly, it can be shown that the other boundary integrals 
namely II, Ib and It corresponding to the left, bottom and 
top boundary, respectively, also satisfy the same inequality as 


that of I R. Therefore, we have 

I|m(x, 7’)||i2(n.) < ||Mo(x)||L2(n,). (62) 

Since T e R+ is arbitrary, one replace T with t in the above 
inequality. This result solution also guarantees the uniqueness 
of the solution. 

For the high-frequency approximation, we have the follow¬ 
ing equivalent formulation on Q, ; Putting F = 5Q, and F^, the 
set of four corner points, define the Laplace-Beltrami operator 
as Ar so that 

idtu + All - 0, (x, t) e Qi x R+, 

m(x, 0) = Mo(x) e suppMo(x) c Q,-, 

dnu + = 0, x e F \ Fc, 

+ diMvj + = 0, {xi,X2) 6 F,- Pi F j, 

where i + j, and i,j e {r,t,[,b}. Taking equation ( l59b and 
following the standard approach, we have 


r (uNu) 

■ dq 


Jr 

j 

r udl^^u\dx\ - j 

r u*Ard, 

J 

r 2 J 

r\rc 

j 

r u*dyu\dx\ + j 

r (druYidrdJ^^^u] 

J 

r 2 J 

r\rc 


(64) 


3 

—i 

4 


Taking the real part and plugging the result back into equation (l59l) yields 

I 


I|m(x, 7’)||^2(n.) - l|Mo(x)||^2(n,,) 


= 2% f dt\-e'’'l^ f u^Sy^dxl + f (druy(drd;^^\)ldxl 

Jo L Jr 2 Jr\rc 

= 2 r d^y f dd-\^\^iy\\^\-^iyA\^Mu{x,m\ 

-l|f<o) p Jr ' ^ > 


(65) 


where q is taken as the tangential variable to the boundary so that dr = d^. In the region defined by |^| > |^P/2, we obtain the 
condition 

I|m(x, T)||/,2(n,) < ||Mo(x)||/.2(n,), (66) 


which guarantee stability and uniqueness of the solution on Schrodinger equation given by 
the bounded domain. 


III. GENERAL SCHRODINGER EQUATION 

The linear Schrodinger equation with time-dependent po¬ 
tential and the nonlinear case are treated in the same fash¬ 
ion. Let us consider the IVP corresponding to the general 


idtU + Au + (f>(x, t, \u'^)u - 0, (x, f) e X R+, 
m(x, 0) = Mo(x), X 6 R^. 


( 67 ) 


The potential function, 0 is assumed to be real valued. 
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FIG. 5. Parametrization of the boundary F = 30, 


A. Convex domain with smooth boundary 


The computational domain is taken to be a convex set 
fi, c with a smooth boundary T = 30,. Let be the 
parametrization of the curve T where s is the length along the 
curve. Introducing the tangent vector e^, also a function of s 
only, and taking into account the convexity of the domain we 
have 


dx — dre„ + (1 + rK)dstr, 


( 68 ) 


where k is the curvature of the boundary T given by 


k(s) = 


d&T 

ds 


|det(y',y')| 


(69) 


The initial data is assumed to be compactly supported in the 
computational domain O,. Carrying out the gauge transfor¬ 
mation u — ve'®, the evolution operator L = idt + h + (pm the 
curvilinear system is given by 


L{r, s, t, dr, ds, d,) - idt + df +Adr + h ^ds + Bds + C, (70) 


where A = (2/(1), H- /j-'x), B = (2ih-^(t>s + h-^dsh-^) and C = 
/A<1) - (Vch)^. The pseudo-differential approach allows us to 
construct various order ABCs given by 
ABC la: 


d„u + - /3^)'^^e"'®M = 0; (71) 

ABC2a: 

d„u + 

+ ^Ku + - id])-^'^ds(e-''^u) 

+ i//ce'®(3, - /3j)“'3^(e“'®M) = 0. (72) 


where 


(d, - id]fi^f(s, f) = (d, - id]){d, - idl)-^i^f(s, t), (73) 


assuming f(s, 0) = 0. 

Remark 3. For the field u(x, t), the action of operators of the 
form (d, — id^ff,a — 1/2, —1/2, — 1,..., can be easily com¬ 
puted by observing that the field can be given a periodic ex¬ 
tension in terms of the parametrization variable, s, so that the 
method discussed in Remark\J\can be employed. The smooth¬ 
ness of the boundary is therefore a necessary condition for 
this method to be applicable. From the previous sections, we 
know that the operator represented by (dt — id^)^"’^^, where m 
is a positive integer, is defined as 

(d,-id])-'"i^f(s,t) 

(74) 


Let 7^sf(s,f) — f(^,f), then the operation (dt — id^) in¬ 
volves the inverse Fourier transform of /(/", t')e~‘^ ^ with 

respect to C,. Therefore, if f(s, t) is of the Schwartz class 
(with respect to s e M.) then so is f(^,f)e~'^^‘~‘ \ so that 
(dt — /3^)“"'^^/(s, f) will also be of the Schwartz class (with 
respect to s e R). Further, it is straightforward to show that 
(dt — id^)^"'^^f(s, f) is continuous at t - 0. 

Next, our aim is to define the operator (dt - /3^)“™^^ for 
tempered distributions f(s, t) such that supptf C [0, oo). To 
this end, let us observe that the transpose '(dt — id^)~'"^^ is 
given by 


'(dt-id^,)-’''/^g(s,t) 

(75) 

where denotes the Weyl fractional integral. The do¬ 

main of definition of this operator is evident. Now, for 
any g(s,t) 6 ,5^(R X R) with supp,g C [0, (x)), the opera¬ 
tion (dt — id^)~’"^^f(s, t) for distributions f(s, t) e c5^'(R x 
R) with suppt f c [0, oo) can be defined by introducing 
a Schwartz class function gi(s,t) such that it agrees with 
'(dt - id])~’”^^g(s, t) for t e [0, oo) so that 

{g(s,t),(dt-id])-'''^^f(s,t)) 

^(:(d,-id])-'''l\(s,t),f(s,t)) (76) 

= {g\(s,t),f(s,t)). 

Therefore, (dt - /3^)“"‘^^/(s, t) is also tempered such that 
suppj (dt — id^)^'"^^f(s, t) C [0, oo). Note that any function 
of the tangential variable, s, can be extended periodically on 
the whole line. Thus periodic functions being a tempered dis¬ 
tribution can be easily included in the domain of definition of 
(dt — id^ff'"^^. It is interesting to note that by applying the 
Leibniz formula for the fractional integrals one can obtain lo¬ 
cal approximation of the ABCs with respect to s: 
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The second family of various order ABCs are obtained un¬ 
der high-frequency assumption with respect to the temporal 
frequency. These are given by 
ABCifc! 

dnu + = 0. (78) 

ABCa/,: 


d„u + e + —ku 


1 


2^ 


,m/4 


a2 

4 ^ 


= 0. (79) 


ABCa*: 


Q , —IkIA /Oi:il/2 —/O . ^ 

Onii + e ' e oJ e u -\—ku 

t 2 


2^ 


,i7r/4 


^ 'i2 

4 




+ — 
2 


/ r A-^ H- d, 


- Y + ds{Kds) 


e®5r'e-®M = 0. (80) 


7. Stability and uniqueness 


In order to study the stability property of the solution of 
the IVP defined by (l67l) with boundary condition ABCia and 
ABC 2 a, respectively, we start with the relation: 


I|m(x, 7’)||^2(n,)-||Mo(x)||^2(n,) 


= 2)R 




(uiVuydg. (81) 


Realization of the boundary operators require an appropriate 
Fourier representation with respect to the tangential variable, 
s. To this end, one may either employ the Fourier series rep¬ 
resentation by extending the field periodically for all i 6 R or 
Fourier transform representation by extending the field to all 
i 6 R by setting it zero outside F. The result for the first oder 
ABCs, ABCifl, can be obtained by observing that the bound¬ 
ary integral 

dt £ (Me-®) - id2(ue-‘‘^)ds, (82) 


satisfies %T i < 0 so that 

I|m(x, 7’)lli2(n,) < l|Mo(x)||^2(n,)- (83) 


For the second order ABCs, ABC 2 fl, the energy estimate can¬ 
not be obtained for a general potential function. A special case 
of interest is when d^ip - 0 so that ABC 2 a is given by 


dnU + e '®m + -m 


+ -iKe'^^Ot - id^) *5^(e '®m) = 0. (84) 


Define the boundary integrals 

-T 


J2 = 


K(s)\ufds, 


I 3 


'iIH 

- ^ C dt f/(:(s)(ue^‘'^)(d, - id^y^d^(ue~"^)ds. 
2 Jo Jr 


(85) 


It follows that %T 2 = 0. Define the symbol crp — 
(1/2)a'(s)(7^ + so that the symbol for the adjoint 

operator Cpt works out to be 


CT pt 




( 86 ) 


Let 2g = P -I- so that 


2ac 


i 


2^ + e 
1 




km 

■2 


(87) 




-t- . . ., 


and putting if/ - ue '®, we have 


X (88) 

From the asymptotic expansion, to the zeroth order, we have 
%!■}, =» 0 (for constant curvature boundary, this result is ex¬ 
act). The energy estimate can thus be approximately estab¬ 
lished for the special case when dstj) - 0. 

For the ABCs obtained in the high-frequency approxima¬ 
tion, it is evident that the first order ABC, ABCi^, satisfies the 
energy estimate. For the higher-order ABCs, one requires an 
additional conditions, dstp = 0 in order to obtain the energy es¬ 
timate. In ABC 2 Z,, the term with the factor a'/ 2 can be ignored 
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as it leads to a purely imaginary quantity. Putting i/f = ue'®, 
consider the following boundary integrals 

Ji= - r dt f 

Jo Jr 2 

-^2 =f dt ^^'‘{s,t)i^&,^^^tf/{sj)ds. 

8 Jo Jr 

(89) 

The hrst integral can be dealt with in the manner done be¬ 
fore. The second integral can be dealt with in the manner 
described in the appendix. Dehning trp - so 

that cTpt = Putting 2Q - P + we have 

ctq - -/c^ cos[; 7 r/ 4 -H 7 rsgn(^)/ 4 ]. The boundary integral can 
now be written as 

X W^)]2|^r‘/2^os[;r(l + sgn(^))/4]|.^,[(A(^,f)]l^ 

(90) 


so that %l 2 < 0 . 

For ABCafc, we dehne the following boundary integrals 


12 = 
IA = 


^ f dt f t//*(s,t)[(/<^ + d^/()df ^]i//(s,t)ds, 

8 Jo Jr 

Js =J" dtj"i/r'‘(s,t)[(d„<p)d^Ji/r(s,t)ds. 

Dehning crp - (/c^ -i- 5^/c)(i^)“' so that 

o-Q = + d]K)\^l^ cos[;rsgn(^)/4]. 

From here it follows that J 3 is purely imaginary. Carrying out 
an integration by parts in J 4 with respect to s, we have 


=41^'! 

"■ 4 X 


[dsilJ*(s, t)]Kd^ ^ [ds^{s, t)]ds 


(92) 


KiiO-'\Mdsilj{s,t)rds. 


Therefore, J 4 is purely imaginary. The discussion of the 
last integral is more involved. Dehne the symbol crp - 
-{dn<p)iiO ^ so that the symbol for the adjoint operator cr^t 
works out to be 


Now 


cTpt ~ - ^ ^(5f(9„(^)(9|(-!^) 

teN 

2o-q ~ -{d„4>){i^y^ - ^ ^(5f5„<^)5|(-!^)^' 
teN 


(93) 


( 94 ) 


Therefore, the contribution from this last integral can be ig¬ 
nored on account of the lower order leading term in the above 
asymptotic expansion. Thus the energy estimate 

ll“(X’^)llm.)^ll“o(x)lli 2 (^,). (95) 

can be established for the ABC-family^: ABCj*, j - 1,2,3. 


B. Domains with straight boundary: infinite strip 

The ABCs for the straight boundary can be written by set¬ 
ting the curvature k to zero. 

ABCifl! 


d„u + e'“^^e®(5, - idy^e~''^u = 0. (96) 

ABCia: 


dnu + 

+ - id]r^'^ds{e-‘'^u) = 0. (97) 

Along the direction in which the strip extends to inhnity, we 
impose the periodic boundary condition. Note that an exact 
representation of the boundary operators is possible in this 
case which is based on the observations made in earlier sec¬ 
tions. The energy estimate for the hrst order ABCs is easy to 
obtain; however, for the second order ABCs this result is not 
available. The high-frequency ABCs work out to be 
ABCii,: 

d„u + = 0. (98) 


ABC 2 i: 


dnU + 

1 






¥ 


0. (99) 


ABCs*: 


dnU + e '®m 

( 100 ) 

The energy estimate for the general case is not available for 
such ABCs. 


C. Rectangular domains: comer conditions 

Continuing with the rectangular domain as dehned in Sec¬ 
tioning let us consider the possibility of extending the results 


'* Only the second and third order ABCs require = 0 in order to establish 
the energy estimate. 



12 


obtained in Section IIIB II to the nonlinear/variable potential 
case. We confine our attention to the high-frequency ABCs 
only. We demonstrate the posibility of constructing the corner 
condition at the corner defined by Tr fiF,. Consider the ABCi* 
given by 

- 0 , x 6 ( 101 ) 

dx^u -t- 

- 0, x e V,. (102) 

I 


Let us consider the weak formulation of the IVP defined by 


idfU + Au + (pu = 0, (x, t) 6 Q.i X R+, 

m(x, 0) = Mo(x), X 6 Qj, with suppMo(x) c Q,-. 

Consider the test function i/r(x) e //^(Q,) so that 


(103) 


(idfU + V^M -I- (f)u)il/(fx 


f d' 

Jcii 

- I [id,u - 
Jni 


(Vu) ■ (Vifr) + (pu\(rx + I ij/iyu) ■ dq. 


X' 


(104) 


Let us consider the top and right boundaries. The boundary 
integrals are given by 


Ir + It - ^ il/(dxiU)dx2 + J" i/^(dx2u)dxi 






/ 

Jr^uF/ 

f 

Jr^-uF/ 


' 2 


1 


X 


il/dx2(e‘'^d, '®m) 


X 2 =xt Jr, 


X 


x- 


“^u)dx2+ ( ‘'^u)dxi 


(dx2d)dx2d"^d, ‘^u)dx2 


(105) 


-H i/^5x,(e'®5, '®m) 


- r {dxy)dxid ‘‘^u)dxi 

H=xi Jr, 


From here it is evident that the corner condition must provide 
the value of the following expression at F^ Fi F,; 

5x,(e®57^^^e-®M) -h 5x,(e®57‘^^e-®M) 

= /((fix, +(fix,)e®57'/'g-''®M-/g®57'^2[((fix. +<fix.)e-®M] 
-H e®57'^V®(5x,M -H 5 x,m). (106) 

For the linear case with variable potential, the last term in the 
above equation must be computed from a corner condition. 
However, for the nonlinear case, <fi is supposed to be depen- 

I 


dent on some known field while our final intent is to restore its 
nonlinear dependence. In this light, the above equation can¬ 
not be made free of the derivatives of u at the corner as the 
equation is nonlinear. However, we may derive a condition at 
the corner which can be combined with the ABCs and eventu¬ 
ally be solved by an iterative scheme. Putting u - in the 
evolution equation (I1031 l. we have 

id,il/ + Ail/ + 2/(V(fi) ■ (Vi/r) -H (e-®Ae'®)i/^ = 0. (107) 

Carrying out the operation on both side of the equa¬ 

tion above, we have 


/e®5y V + -n 2e'®57‘^^(/V(fi) ■ (Vi/r) + e® 57 ‘^^(e-®Ae®)iA = 0, 

/e®5‘^V + e®57*^^Ai/r h- 2e®(/V(fi) ■ 57‘^^(ViA) -h (Ae®)5,:‘^V + ■ - ■ = 0, 


F 


which simplifies to 


This equation combined with the ABCs gives the following 
corner condition 

CCi : dx,u + dx2U + = 0, (110) 


ie‘^d]^^e "^u + 


A^e'®5, '®Mj= 0 mod(5, '®m). 

(109) 


or, equivalently. 


CCi : e®<970V®(5x,M -H 5 x,m) -h = 0. (Ill) 


3 _,v 
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Similarly, for ABCa;, one has the following corner condition 
CC2 : dxiU + dxjU + 

( 112 ) 

or, equivalently, 

CC 2 : e‘‘^d;^'^e-'‘^{dx,u + dx,u) + 

= 0 mod 

(113) 


Schrodinger equation and general Schrodinger equation for¬ 
mulated on various types of computational domains to be ex¬ 
pressed in a natural way. In particular, two families of ABCs 
are studied in this paper; first obtained with the standard 
pseudo-differential approach and the second one obtained as 
the high-frequency approximation of the former. For the 
rectangular domains, we have developed various order cor¬ 
ner conditions for the family of ABCs obtained in the high- 
frequency approximation. Each of these families of ABCs 
(along with corner conditions) are also investigated for sta¬ 
bility and uniqueness of the solution of the resulting initial¬ 
boundary value problem. The discussion of the time-discrete 
version of the TBCs/ABCs is beyond the scope of this paper; 
this investigation is in progress and it will be published in a 
forthcoming paper along with the numerical results. 


D. Special case: ^ - (p(t) 


Appendix A: Some properties of pseudo-differential operators 


For time-dependent potentials with no spatial variation, the 
quantity O is purely time-dependent and the ABCs for dif¬ 
ferent domains types can be simplified considerably. It must 
be remarked that the exact form of the transparent boundary 
condition is only obtainable for an infinite-strip with periodic 
boundary condition in the unbounded direction or rectangu¬ 
lar domains. At the appropriate segments of the boundary, it 
reads as 

dnU + - 0. (114) 

For rectangular domains, we consider this problem is more 
detail. At the boundary F^, we have 

dx,u + = 0. (115) 

The auxiliary function, (p(xi,X 2 , t, f) as defined in Section ing 
in the present case is given by 

^X2[<fixi,x2,t,t')] = e“'^2('-'')-'®('')^_^jM(xi,x2,f')], (116) 


Fet us consider the symbol space x R") of M-quasi 

homogeneous symbols S where Y is an open subset of R" 
and M - {jji,fi 2 , ■ ■ is an n-tuple of numbers /r, > 0. Fet 
P be a pseudo-differential operator with the symbol p(y, ^) 6 
S'^(Y X R") so that 


Pu(y) = J p(y, Oe'^'^ii(Od"( 

= JJ p(y,Oe‘^-^^-^'W)dy'd^, 


(Al) 


where u(y) e C[^(T) and d( - d^\d ^2 ■ ■ ■ ,d^n, dy - 
dyidy 2 ..., dy„. The adjoint of this operator, denoted by P\ 
can be defined as 


P\{y) = (27r)-" 



p\y',Oe‘^'^-^'^u(y')dy'd^. 


(A2) 


so that ^(xi,V 2 ,f,f) = u(xi,X 2 ,t)e Putting So(^) = 

=^xi=^x 2 “o(x), we have 


This operator belongs to a more general class of pseudo¬ 
differential operators defined by 


ip{xi,X 2 ,t,t') = f ^ ^ 

(27ry Jr2 

From here is follows that the IVPs satisfied by this auxiliary 
field is not different from that described in Section ing 


Pu(y) = 



p{y,y',Oe'^'^-^'^u{y')dy'd^, 


(A3) 


where the symbol p(y,y', 0 e C°°(T x T x R") lies in *S^(T x 
Y X R"). The adjoint of this operator is defined by the symbol 

p^(y,y',O^P*(y',y,0- (A4) 


IV. CONCLUSION 

In this paper we have discussed the formulation of the op¬ 
erator (dt - /Ar)“, a - 1/2, -1/2, -1,..., in terms of the frac¬ 
tional operators. This allows the TBCs/ABCs for the free 


Note that the primed variable is to be integrated over. These 
new operators do not have unique symbols but they do admit 
of a representation in terms of the former kind of operators as 
P - OP(crR(y, /■)) + R where crp(y, {) e S'^(Y x R") and R is 
an operator with kernel Kj{(y,y') e C“(T x Y) such that 
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a-p(y, 0 


Z 

o'eN", |Q'|<iV-l 


1 


5 “[ 5 Xy,y, 0 ]/= 


6 X y X R"), 


(A5) 


where N > 0 and fi is smallest element of M. The relation¬ 
ship in equation (IA51 l defines the asymptotic expansion of the 
symbol and we write 

o-p(y,0 ~ Yj (A 6 ) 

aeN" 


If crgiy, is independent of y, then 

%{u\Pu) = (2nr f d^imOl^crgiO. (A12) 

The sign of this expression can solely be decided by the sign 
of o-g(0- If crg(y, 0 = (^(yi)cr(^) where y = (y 2 ,y 3 , ■ ■ ■ ,y«) e 
R"-i and ^ = (^ 2 ,. ^«) 6 1 ^"“' then 


This determines <Tp(y, only up to a smoothing operator. The 
formula expressing crpf(y, in terms of cr* (y', can be ob¬ 
tained by writing 


%(ulFu} = 



dyi dl\^yu(y \, y)|^<^(y i)o-(^). 


(A13) 


p*(y',y, - a-piy', 0 e S^-(Y xYx R"),, (A7) 


and using the asymptotic expansion 


<rpt(y, (A8) 


cteN” 


where R being a smoothing operator does not show up in the 
asymptotic expansion. The adjoint can be used to write the 
Fourier transform of Pu(y) by noticing that (Pt)t = p^ so that 

PM(y) = (2;rr J(J cr;,(y',Oe-‘^'^''u(y')dy')e‘^'^d{. (A9) 

Consider the inner product defined by {m|v) = f u*vdy. An 
expression of the form '’K{u\Pu) arises in establishing the sta¬ 
bility of the ABCs. This can be computed by observing 


2%{u\Pu) — (mIFm) + {u\P^u). 


(AlO) 


On account of the ambiguity in the knowledge of the exact 
symbol, the conclusion remains valid only up to an infinitely 
smoothing operator (except for the special cases where the 
symbol is of principle type). 

• Fractional operators with symbol p(t,^) = 

(where t e R with covariable ^ 6 R \ {0)); 




= i7ii^rcos 


—sgn(^)H-arg ?7 


(A14) 


• Operators with symbol p(x, t, = rjii ^+(where 
(x, f) e R^ with covariables (^,^) 6 R^ \ {(YO 6 • 

f + ^2=0)): 


~ + irr + - iO"] 


= |77l|^H-^Ycos 


na ^ 9 ^ 

—sgn(f H-^ ) + arg ;7 


(A15) 


• Let cTp = rj^(x)(i^Y with real valued function (f>(x) > 
0, Vx e R where the meaning of the variables is same as 
that of the last example. Then 


Define 2Q = (P + P^) and using Plancheral’s theorem, we 
have 


= (All, 


crgix,t,Y^) ~ -mxmr + p<p(x)(-im 




/.-X 

—sgn(^)-Harg ;7 


(A16) 


so that 




I// 


d^dx\^t[u(x, t)](x, ^)r 0 (x)|^r cos 


na 

—sgn(f) H- argT/ 


(A17) 


In all the above cases, when a < 1 and arg rj - njA, it is easy exactly in the main body of the paper for the aforementioned 
to show that the sign of erg remains fixed. This case is treated operators without resorting to the properties of the pseudo- 
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differential operators (i.e. the ambiguity resulting from the lack of knowledge of the exact symbol is circumvented). 
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